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Abstract 

The anisotropic fluid equations governing a frictionless obliquely-flowing plasma 
around an essentially arbitrarily shaped three-dimensional ion-absorbing object in a 
strong magnetic field are solved analytically in the quasi-neutral drift-approximation, 
neglecting parallel temperature gradients. The effects of transverse displacements 
traversing the magnetic presheath are also quantified. It is shown that the paral- 
lel collection flux density dependence upon external Mach- number is nooCsexp[— 1 — 
(M||oo — M_L cot 9)] where 6 is the angle (in the plane of field and drift velocity) of the 
object-surface to the magnetic-field and M||j^ is the external parallel flow. The perpen- 
dicular drift, M_L, appearing here consists of the external E A B drift plus a weighted 
sum of the ion and electron diamagnetic drifts that depends upon the total angle of 
the surface to the magnetic field. It is that somewhat counter-intuitive combination 
that an oblique (transverse) Mach probe experiment measures. 

1 Introduction 

Ion collection by solid objects immersed in a plasma is a problem of perennial interest in 
plasma physics. It provides the basis for the measurement of plasma parameters by electric 
(Langmuir) probes[l| as well as the charging of dust[2| and spacecraft!^. The present work 
addresses the situation where the ion Larmor radius (in the background magnetic field B) 
is much smaller than the object, so that perpendicular plasma fiow is strongly constrained. 

This problem has important similarities to the solution of the fiow of plasma to a plane 
aligned obliquely to the field, the most obvious example being a tokamak divertor plate. 
That problem can be formulated[4] as one-dimensional, taking the coordinates in the plane 
as being ignorable. However, it is well established that no spatially-varying solution of the 
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quasi-neutral plasma equations in one dimension is possible without additional sources of 
particles (e.g. through ionization) or momentum (e.g. from collisions). Recent studies of this 
problem of the one-dimensional magnetized plasma and oblique presheath (e.g. 0, [sl) have 
mostly focussed on collisions as the mechanism allowing the acceleration of the plasma into 
the magnetic presheath. For localized probes, however, if one conceptualizes the problem 
as being dominated by the one-dimensional dynamics along the field, the cross-field flux 
divergence is the most natural effective source to permit parallel gradients. 

Prior theoretical probe studies have focussed on situations where the cross- field magnetized- 
ion flux can be described (somewhat phenomenologically) as diffusive. The full numerical 
solutions for this formulation 0, Isj yield the dependence of the collected ion flux density 
on the plasma density and temperature, and the parallel (to B) Mach-number. That pro- 
vides the theoretical calibration factor for a (parallel) Mach-probe (a probe with electrodes 
facing parallel and anti-parallel to the field), when the perpendicular drift velocity is ignor- 
able. This "calibration" proves to be in good agreement with independent measurements 
and calculations [9j] and has been widely adopted for experimental interpretation. 

The approximate one-dimensional diffusive treatment has been generalized lOl] to include 
an additional perpendicular plasma drift velocity, accounting for the boundary-condition 
modification [111] that the transverse drift causes. Measuring the dependence of the ion col- 
lection current- density on orientation of oblique probe faces then allows one to deduce the 
perpendicular as well as the parallel external drift velocity. The generalized solution can 
be shown|l| to be a simple Galilean transformation of the solution for zero transverse drift, 
which incidentally reminds us of the elementary physical equivalence of E A B drift past a 
fixed object and motion of the object through a stationary plasma. The equivalence also in- 
dicates, though, that the generalized diffusion solution is rigorously vahd only for an oblique 
surface of effectively infinite dimension in the transverse drift direction (so that the Galilean 
equivalence in this direction is valid), but finite in the direction perpendicular to both flow 
and magnetic field (so that diffusion in this perpendicular direction dominates the cross-field 
divergence). Practical Mach-probes generally do not have this configuration. They are more 
often multi-faceted 'Gundestrup' types 121. Il3l. Il4l]. where many short adjacent collectors are 
used with different orientations. So it is not obvious that the generalized diffusive solution 
applies. 

In fact, when there is substantial pre-existing cross-field drift of the ions, it is perhaps 
physically more reasonable to regard that drift as the dominant cross-field transport mech- 
anism, and to ignore diffusion. Gunn 15, loj has explored this problem, with a uniform 
impressed cross-field drift, in two dimensions with his particle-in-cell code. This drift physics 
is appropriate to many space and astrophysical problems too; for example to the interaction 
of Jupiter's satellites with its magnetosphere. It is the purpose of the present work to de- 
rive a general analytic (3-dimensional) solution to this purely advective problem, with fully 
self-consistent drift velocity. 

First, to introduce the solution by characteristics, we recall a recent complete exact 
solution to this problem 16| for an arbitrary shaped probe under the model ansatz that the 
perpendicular drift velocity is uniform. This is a generalization of an earlier self-similar 
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solution 171] mathematically equivalent to a one- dimensional free expansion into a vacuum. 
The very simple general analytic result obtained for the ion flux is gratifyingly close to the 
diffusive-plasma result, and hence to the PIC calculations of Gunn (which include full ion 
distribution-function parallel gradients). The solution demonstrates that provided the probe 
is convex, the flux is not affected (for negligible Larmor radius) by spatial derivatives of the 
surface angle in the drift-direction. The uniform drift ansatz is justified by inspection only 
when the probe is two-dimensional; so that the coordinate perpendicular to the field and 
drift is ignorable and the probe-perturbation of the plasma does not introduce additional 
drifts except along the ignorable coordinate. Again, practical probes are generally not well- 
approximated as two-dimensional, so the question remains as to whether that uniform-drift- 
velocity solution applies in practice. 

The following remarkable result is rigorously demonstrated in section |H The ion flux to 
the probe surfaces derived for uniform-drift-velocity does apply even when the full spatially- 
varying self-consistent drift velocity, including the perturbation from an arbitrarily-shaped 
three-dimensional probe, is accounted for. When the external drift arises purely from electric 
field, one can obtain the full self-consistent spatial dependence of the density and velocity 
throughout the perturbed plasma region, using an elementary geometric algorithm. Some 
examples are given. 

Furthermore, external diamagnetic drifts can also be included, again for arbitrary-shaped 
three-dimensional probes. They make important but counterintuitive contributions to the 
observed ion current density. In addition to the effects that arise in the plasma, it is essen- 
tial to account for transverse displacements that arise in the magnetic presheath; they are 
calculated in section [5l Such local drifts in the magnetic presheath have previously been 
identified [l8|, [l9| as important contributors to oblique boundary conditions. Based on these 
prior considerations, it has sometimes been asserted that diamagnetic drifts do not trans- 
port mass to a surface. The present work provides a more general solution of the magnetic 
presheath displacement effects, dispensing with small-angle approximations, and clarifies 
rigorously the extent to which diamagnetic drifts actually do contribute to the observed 
flux. 

The final result is that a transverse Mach-probe measures effectively the sum of the 
external E A B drift and a combination of the ion and electron diamagnetic drifts. At small 
angles between the field and the collector, the dominant diamagnetic term is the electron 
diamagnetic drift, which of course is generally in the opposite direction to the ion diamagnetic 
drift. 

The presheath displacements can give rise to bias in Mach probe measurements. Its 
relative magnitude is of order the ratio of Larmor radius to electrode size. The effects of 
orthogonal displacements in the plasma region are also calculated rigorously. They mod- 
ify the expression for the flux in ways that are usually of little importance for practical 
measurements. 
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2 Formulation 



We analyse the dynamics of the ion-fluid through the steady-state continuity and momentum 
equations 

V.(nv) = (1) 
mn(v.V)v = —nZeV4> — Vp + nZe(v A B) , (2) 

where m, Z, n, p, Ti, v are the ion mass, charge-number, density, pressure, temperature, 
and velocity, and is the potential. We split the momentum equation into the components 
parallel (y) and perpendicular (x) to the (assumed uniform) magnetic field, and take the 
cross-product with B of the perpendicular part to obtain the form 



(1 \ TTl 



B 

A^J. (3) 



We can immediately identify the first two terms in this expression as the E A B and dia- 
magnetic drifts. The last term can be considered to be the polarization drift, which we will 

regard as ignorable. The approximation of omitting the polarization drift requires the Lar- 
mor radius to be small c.f. the perpendicular scale-length, generally the probe dimensions. 
It can be shown by a posteriori calculation that the polarization drift is smaller than the 
imposed perpendicular drift by a factor that is second-order in the Larmor radius. Ignoring 
the polarization drift term is the meaning here of the expression "drift-approximation" . By 
taking B to be uniform we have of course eliminated the grad-B and curvature drifts. We 
adopt the simplest possible fluid closure scheme, that the ion temperature, Tj, is invariant, so 
that the pressure is simply proportional to density. Together with dropping the polarization 
drift, this makes Vj_ divergenceless: 

Under these approximations the continuity and parallel-momentum equations become 

(v.V)lnn + V||W|i = (5) 
(v.Vh + ^V||(0+|^lnn) =0, (6) 

while the perpendicular-momentum conservation is expressed by the drift expression. 

The potential is eliminated from these equations by accounting for the self-consistent 
solution of the electric field arising from the ion and electron densities. The electron density 
response along the magnetic field is, as usual, taken to involve rapid equilibration; so that 
the electron pressure gradient is balanced by electric field: 

V||0=(re/e)V||lnne, (7) 

where the parallel gradient of the electron temperature, VyTg, is taken to be zero and the 
electron density is rig. Assuming that the Debye length is much smaller than the probe, we 
will treat the plasma as quasi-neutral, so that rig — Zn. 
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Using the notation = {ZTf, + Ti)/m and M = v/c^, the ion equations then take the 
normahzed form: 



M.Vlnn + V|jM|| = (8) 

M.VM|| + V||lnn = 0, (9) 

which can be rearranged by adding and subtracting to show exphcitly the "characteristics" 

(M.V + V||)(lnn + M||) = (10) 

(M.V- V||)(lnn-M||) = 0. (11) 

Thus the quantities (Inn ± My) are constant along their corresponding characteristics = 
(M ± 'Q/B)ds. And we can fully solve the problem by analysis of the characteristics. 



3 Uniform perpendicular-velocity ansatz 

First we review the solution under the condition that the perpendicular velocity is simply 
a constant, M^, independent of space. This ansatz is clearly justified if the coordinate 
perpendicular to B and is ignorable. See ref 16|| for additional details and explanation 
of the following derivation. We choose axes such that B is aligned along x and = M^y 
along y. The requirements expressed by the characteristics (ITOl [TTl) are that both 



(Inn + M||) = const along dx = dy{M\\ + l)/Mh, (12) 

and 

(lnn-M||) = const along dx = dy{Mn - 1)/Mh. (13) 

These will be referred to respectively as the positive and negative characteristics. 

For definiteness, we now consider plasma that is on the higher-x side (to the right) of the 
object. Figure [T] will be used for illustration. For any point in the plasma, a positive and a 
negative characteristic pass through it. If both characteristics originate in the unperturbed 
plasma at y ^ — oo, and do not enclose the object, then values at the point satisfy both 
Inn + M|| = Innoo + Moo and Inn — My = Innoo — M\\oo (where oo indicates values in 
the unperturbed region). These simultaneous equations have only one solution: n = n^o, 
M|| = M||oo, showing that the point is in the unperturbed region, for example: Pq. The 
characteristics for such points are straight lines with slopes M/i/(M||oo ± 1). 

The most important case is when just one of the characteristics originates not at y = 
—oo, but on the surface of the object (e.g. Pi). The positive characteristic is always to 
the left of the negative characteristic, when tracing backward from a common point. So 
the characteristic that originates on the object is the positive one. On that characteristic, 
(Inn + M||) is constant, but not equal to the unperturbed value. Each point along the 
positive characteristic satisfies also Inn — M\\ = Innoo — ^||oo because there are negative 
characteristics from infinity to each point. The only way to satisfy these two requirements 
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Figure 1: (Color online). Construction of the solutions at different points. Pq is in the 
unperturbed region. For Pi the characteristics are shown. P3 is in a concave region and 
so its positive characteristic is not tangent at P3. A value 6 > 7i/2 such as for P4 is not 
problematic. 

is that, along the positive characteristic, My = const and n = const. If My = const, then 
the slope of the characteristic, M/i/(M|| + 1) is constant. It is a straight line. 

The line's slope is determined by the absorbing boundary condition at the plasma edge. 
That condition requires |16i] My to be as negative as possible consistent with the overall 
solution, which requires the greatest possible slope-angle 9 = arctan[M/j/(M|| + 1)] (even 
perhaps such that 9 > 7r/2). The characteristic must therefore always be tangential to 
the object boundary where it intersects it. Thus, all positive characteristics that originate 
on the boundary do so as tangents, and for any point in the perturbed plasma region the 
positive characteristic is that straight line passing through the point which has greatest 9 and 
originates as a tangent on the object. Once that line is determined geometrically, its slope 
determines My and hence n at all points along it. If the steepest tangent angle is less than 
9m = arct an [Mft/( My 00 + 1)], then the positive characteristic does not intersect the object, 
but has slope 9 = 9m', and the point is in unperturbed plasma. The entire solution for the 
plasma in the perturbed neighborhood of an arbitrary-shaped object is thus 

n = rioo exp(My - Myoo), My = Mh cot ^ - 1. (14) 

The solution f|T^ provides an extremely simple formula for the ion flux to a surface not 
affected by concavity. Adding the perpendicular and parallel components, the total flux 
density along the outward normal within the {x, ?/)-plane, i.e. in the direction (— sin^, cos^), 
is nCs{Mh cos 9 — My sin ^) = ncg sin 9. Written as flux per unit area perpendicular to the 
magnetic field, this is 

Fy =nooC,exp[-l - (My^o - M,, cot ^)]. (15) 
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First, this form indicates, importantly, that for points not in a concave region of the 
object, the collected flux depends only on the angle of the surface there, and not on the 
shape of the object at smaller y. Second, the exponential dependence upon M\\oq is within 
10% of the dependence, exp[— 1 — l.l(M||oo — Mh cot 9)], that fits the diffusive solution[lol. 
Third, consideration of the characteristics shows unambiguously that leading faces, for which 
6 < 6m receive simply the unperturbed flux [Fy = nooCs{Mh cot ^— Myoo)], while trailing faces, 
even those for which Myoo — cot 6 > 1, are governed by the formula f[T^ . The boundary 
condition at the magnetic presheath edge is just the same as what is sometimes called the 
"magnetized Bohm condition" but arises here naturally from the analysis of the quasi-neutral 
equations. 

In a concave region of the object, e.g. at P3, the surface angle (local tangent) 9s, is 
smaller than the characteristic's angle 9, and the distinction must be retained. This leads to 
an enhanced ion flux, equal to equation fllSp times the extra factor {Mh cot 9s — Mh cot 9 + 1). 

4 Accounting for self-consistent drifts 

The presence of the probe perturbs the plasma potential and density, for example as cal- 
culated in the model case of the previous section. Most probes are of limited extent in the 
direction {z) perpendicular to the plane containing B and M^oo, and indeed may have a 
non-zero 2;-component of the normal to the collecting surface. The plasma perturbations 
therefore give rise to spatially-varying ion drifts that frequently break the (z-translational) 
symmetry assumption used to justify the homogeneous-M^ ansatz. So we must now return 
to the full equations ( ITOl [TTj) accounting for the complete, self-consistent, spatially-varying, 
perpendicular velocity. 

The interesting case is of points for which (only) one of the characteristics starts not at 
infinity but on the probe itself. As before, for definiteness, but without loss of generality, we 
will take that to be the positive characteristic. At this point then (Inn-t-My) 7^ (lnnoo+M||oo); 
nevertheless, because of the negative characteristic, it is still true that (Inn — My) = (Innoo — 
^Ijoo) = const. This shows that in a region whose points have any characteristic starting 
at infinity, the nature of the solution is of the form My = M\\{n), and especially VMy 
is parallel to Vn. Consequently if there is a self-consistent combination of density and 
velocity fields {n, My, = M^q} that satisfies the advection equations f[TU| [TTl) (and 
the drift equation, H]), any perpendicular vector field, M_|_i, that satisfies M_|_i.Vn = 0, also 
satisfies M_li. VMy = 0. We can therefore subtract any such M_li from M_l without affecting 
the characteristic equations (fTOl [TTl) . In other words, the combination {n, My, M_l = 
(M_Lo — M_Li)} also satisfies the characteristic equations (though not the drift equation, 
H]). Moreover, the subtraction of M_|_i leaves the boundary condition at the plasma edge 
invariant provided probe curvature is small compared with 1/ps- The invariance follows 
from the boundary condition being that the positive characteristic be tangent to the surface. 
If the surface is expressed by the equation s{x, y, z) = 0, then tangency is cix.Vs = 0, 
which along the characteristic is (My + l)Vys -|- M_l.V_ls = 0. In so far as Vs and Vn are 
parallel (i.e. the density is invariant in the tangential directions) subtracting M_|_i leaves this 
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condition unchanged. But n is indeed invariant along the surface when given by eq 
provided that the surface angle {Og) gradient (i.e. curvature) can be ignored. 

The significance of these observations is profound. It means that although the actual 
perpendicular velocity may be very complicated, and include drifts arising from the 
self-consistent potential- and density-gradients caused by the presence of the probe, we do 
not have to solve the associated complicated characteristics. Instead, we can subtract from 
any drift that satisfies M^i.Vn = 0, and solve the resulting simpler characteristics. The 
resulting solution for n and Mn is correct then also for the full drift expression. 



4.1 Pure E A B external drift 

Consider first the case when the unperturbed plasma has uniform density (rioo = const) 
as well as temperature, but has a uniform impressed perpendicular electric field in the z- 
direction perpendicular to the field , (B/i? = x), giving rise to an EAB drift. In other words, 
we have a non- uniform potential, (j)oo{z), such that V^oo = —zE; so that the homogenous 
drift is = yE/B = MhCs. 

In the presence of the probe, the electron parallel momentum (force-balance) equation 
can be integrated along the field, from infinity to any position to give the perturbed potential 

- 4>oo{z) = (Te/e) ln(n/noo) . (16) 

We substitute this into the drift expression to get 

v± = -V U^ + -^\n{n/n^) + j-\nn] A— (17) 

= -(v0oo + ^Vlnn^ A^. (18) 

Or, in normalized form: 

M± = Mh- Ps V In n A B/B , (19) 

where ps is the ion Larmor radius at the sound speed. The perpendicular drift veloc- 
ity thus consists of a uniform term M/^, equal to the external drift, plus a term M_li 
(= — psVlnn A B/B) arising from local gradients of density (and associated potential), 
which satisfies M^i.Vn = 0. Our approach is therefore to solve along characteristics defined 
not by the complicated full drift velocity M^, but using the uniform external drift velocity 
Nlfi which arises from subtracting off M^i, as we have shown we are permitted to do because 
M^i.Vn = 0. But such an approach, of using a uniform impressed perpendicular drift, is 
precisely the ansatz solved in the previous section, albeit without this detailed justification. 
Therefore, the solution obtained there applies without modification. What is modified is 
that the condition of translational invariance in the orthogonal (z-) direction is removed. In 
other words, the restriction that the probe be two-dimensional, which was invoked previously 
to justify the neglect of the density-gradient-induced drifts is proven here to be unnecessary. 
The results of that section, the dependence of n on Mm , the spatial variation of Mm embodied 
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in the equation tan6' = Mh/{M\\ + 1), and the surface flux expression, apply to any shape 
of three-dimensional probe, provided only that the condition of convexity (that the surface 
is not reentrant) is satisfied. This convexity condition must be applied along the original 
characteristic (before subtracting M^oo) so it depends on the full 3-D shape of the probe 
(and the z-drifts). 

The fiux-density of ions from the plasma to the probe surface in this case is given by 
the vector sum of M/^ and My because the additional term M^i is always locally tangential 
to the surface, even if the surface- normal has a non-zero z-component. Therefore eq f[T^ is 
valid. 




Figure 2: (Color online). Contours of My + 1 = 1 + M\\^ + hi{n/nao) at intervals of 0.1, 
in planes of constant-x near a sphere of unit radius in a plasma with external drift velocity 



Ml, 



0.2, Mh = 0.24. 



As an illustration of this complete solution of the problem. Fig [2] shows a representation 
of the 3-Dimensional variation of My and equivalently In n by contours of the quantity My + 1 
drawn in perpendicular planes at various distances from a spherical object. Regions empty 
of contours to lower-y (left) of the contours shown have uniform unperturbed plasma. To 
the right of the contours is the wake region where the equations are not valid. 

A second illustration is in Fig [31 which shows the contours for a pyramid shaped probe 



similar to what is used in Alcator C-Mod experiments (2 II] . A different drift velocity is 
illustrated. 
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Figure 3: (Color online). Contours of M\\ + 1 = 1 + M\\oo + ln(n/?T,oo) at intervals of 0.1, in 
planes of constant-x near a pyramidal Mach probe in a plasma with external drift velocity 
M||oo = -0.1, Mh = 0.15. 

4.2 Inclusion of external Vn diamagnetic drift 

The presence of a diamagnetic drift arising from external density gradient leads to several 
complicating factors. We consider an unperturbed density, noo, that in this case is not 
uniform but is a function also of perpendicular position. The usual case has V±0oo and V ±noo 
approximately parallel to each other, because both are perpendicular to the flux surface in 
a confined plasma. This is actually essential for the consistency of the unperturbed state, 
if it has no parallel gradients. Then the total drift velocity does not satisfy V.(?7.ooVoo) = 
unless B.(V0oo A Vrioo) = so that the drifts arising from the cf) and n gradients are parallel 
to each other. We stick to this case here, and consistently take all external gradients to be 
in the z-direction, but note that effects in tokamak scrape-off-layers where parallel gradients 
are present and give rise to cross-flux-surface drifts are thereby ignored. In other words, we 
are dealing with of negligible parallel gradients in the external plasma. 

We suppose that the external logarithmic gradient of the density is constant: 



where L„ is the (constant) density scale length. Equations (fT6!) and (fT7|) still apply. We 
still can rely upon integration along the negative characteristic to write down a relationship 
between n and Mm, but that relationship is now 




(20) 




(21) 
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where rir 



is the unperturbed value of density at large distance from the probe, but 



at a value, z^, oi z corresponding to tracking backward along the negative characteristic 
from the point of interest. Because function of z when there are diamagnetic drifts, 

the value of In 77,00(^00) depends upon the total z-displacement, 8z = z — Zoq between the 
characteristic's start and the point of interest. Thus it is no longer the case that n = n{M\\). 

Write the relationship between density and My , deduced from the negative characteristic 
integration as 

\n{n/n^{z)) = -M,,^ + M,, - 5z/L^. (22) 

In this expression, 6z is not a constant. We will demonstrate in the following, however, that 
a solution exists in which 6z is a function only of My. So taking 6z = 6z{M\\) we observe 
that ln(n/noo) is also a function only of My. [Here and following we use the notation rioo 
without an argument to denote the unperturbed density at the position of the point, noo{z) 
not noo(-2oo)-] The drift velocity from eq f|T7|) can be written 



V u, 



^Inno 
Ze 



+ — ^Vln n/nc 
ze 



B 



(23) 



in which the first two terms give rise to perpendicular Mach number = Me + M„j, the 
sum of external E f\B and ion diamagnetic drifts. The final term, which we identify as M^i 
is perpendicular to Vlnra/noo and, because Inra/noo is a function of My, perpendicular also 

to VMy. 

Now we write the positive characteristic equation so as to use these expressions: 







d 



fin n + Mil) = — ^ 
^ ^ cat 

d ' 



c.dt 



(inn/: 

(2My - 6z/Ln) = MjLn + 



+ My j + MjLn 
d 



{2\nn/n^ + 5z/Lr.), (24) 



where we use the notation 

d 



c.dt 



= M.V± Vii 



(25) 



for the derivative along the positive or negative characteristics. Because the arguments 
of the derivatives in the last two forms of eq are explicitly functions only of My, we 
can subtract M_li.V from the characteristic derivative without effect. In other words, in 
those expressions, we can interpret the derivative in the alternative version (valid only when 
operating on functions only of Mn) 



d 



c.dt 



± 



(26) 



Now we eliminate ?/-derivatives of In n from the last form of the positive characteristic equa- 
tion using the expression for the z drift velocity: 



M. 



ps 



dlnn dlnn/ric 



dy 



ps 



dy 



(27) 
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to find 



d 



1 ^A4 

Ps , 



1 d 



2 cAt 



5z_ 



(28) 



In the same way, we express the negative characteristic equation in terms of \n.n/n^, rein- 
terpret the derivative as the form fl26l) and then we ehminate x- and y-derivatives of hm/rioo 
using eqs fl25]) and (ITTI). 



d 



c.dt 



Mil 



d 



c^dt 



M|| - 1 

Mil + 1 



-mJ -— + ^ 



d 



c.dt 



Sz 



^ MhM, ^ M, 

Ps Ln 



(29) 



For compactness we define d/dM\\{6z/2Ln) = r; so that 



c.dt 



6z 



Mil.- 



d 6z 
'dMii 2L„ 



c.dt 



Mil. 



(30) 



Eliminate 



Csdt 



M|| between this expression and the Mm form of eq (124|) to obtain 



d 



c^dt 



Sz 
2U 



M^ r 
'2L„ 1-r' 



(31) 



Substituting into eq (129|) and ehminating M^ and the right hand side's My using the 
identity 



d 



bz 



d 



Mil 



d 



'2iL^ I Cgdt 



Mil 



we arrive at the following quadratic equation for r 



1 



Mil 



+ 3 + (M|| - 1) + AL^Mh/ps 

1 — r 



M|| + 1 

We solve this equation, using the simplifying notation 

a = 1/(4 + ALnMh/ps) , u = aMy 

to find 

r = (1/2)[1 + 2u± Vl - 4a + Au^] . 
Then we can integrate to obtain the 2;-displacement 



^ = 2frdM,= ^- 



u + — u\lv? — a 



a - \)\a. [u + d - a + \ 



(32) 
(33) 

(34) 

(35) 

(36) 
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This solution is real only if i — a > 0, which is equivalent to ps/{LnMh) > —1. Since the 
ion diamagnetic drift is M„j = — ps/L„ and Mh is positive by presumption, this requires 
Mni > -Me and Me > 0. 

In Fig m are plotted the solutions for 6z as a function of My for the entire useful range 
of a, i.e. of MhLn/ps = —Mh/Mni- We take that to be constrained by —Me < Mni < Me 
(and Me > 0) which is equivalent to ps/{LnMh) > —1/2. Whether normalized by L„ or by 
Ps/Mfi, these results show that the displacement, 6^, while non-zero, is modest. Moreover 
the curves are quite close to being mirror-symmetric about the line M» = — 1. 
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Figure 4: Solutions for the 2;-displacement normalized to the density scale-length (a) or the 
Larmor radius (b), for the useful range —Me < M^i < Me- It is assumed that M\\oo is zero, 
but the curve shapes are the same regardless of the starting value of Mii. 



The value of 6z together with the negative characteristic integration, eq (12^ provide the 
complete solution for ln(n/?2oo) which is a function only of My. The characteristic equations 
can be considered to be 

(lnn/noo±M||) = -^, (37) 

in which the characteristic derivatives can be taken as eq fl26p . that is, lying in a 2; = const 
plane. The characteristics are now curved, and cannot be constructed directly from the 
geometry, as was possible in the absence of diamagnetic drift. But this does not matter 
for the purposes of obtaining the flux to the probe. The boundary condition at the plasma 
edge is, as before, that provided the boundary is convex, the positive characteristic must be 
tangential to the surface; that is, M\\ -|- 1 = M^ cot 9. 

In practice, the perpendicular velocity is generally deduced from Mach probe measure- 
ments effectively by comparing values of the ion current density for faces having equal and 
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Figure 5: (Color online). Schematic illustration of shape of positive characterstics in the x-y 
plane. In (a) L„ is positive corresponding to diamagnetic drift opposing E A B drift, and 
the displacement 6i being down the density gradient. In (b) Ln is negative. 

opposite values of cos 6'; that is, having the same angle to the magnetic field, but pointing 
upstream or downstream with respect to the perpendicular flow. Although 6z changes the 
value of In n and thus affects the ion flux density, if 6z is of exactly even parity in My + 1 
and hence in cot 6, it will contribute nothing to the ratio of the ion fluxes on which the 
measurement is based. For L„ —>■ oo (zero diamagnetic drift) or M^L^ = (zero net drift) 
6z is indeed exactly of even parity. The maximum value of the odd-parity part of dz/Ln is 
approximately 0.02 for |M|| + 1| < 0.6 for \Mni/ME\ < 1 (and rapidly decreasing as \Mni/ME\ 
decreases). Thus the contribution of the Sz term to Mach probe velocity measurement is 
typically less than 2%, which is for practical purposes negligible. The magnitude of 6z is 
sufficiently great, however, that it would be inadvisable to attempt to deduce the perpendic- 
ular flow without taking advantage of its approximate parity. In other words, a Mach-probe 
measurement really must be based on electrodes with equal and opposite values of cot 9, and 
not, for example, on comparing positive cot^ with cot 6 = 0. 

It is, in effect, the 6z variation that causes the positive characteristics to have curvature. 
This curvature has magnitude approximately 1/Ln- When considering whether or not a 
probe surface is convex, this curvature has to be accounted for. Schematic representations 
of the positive characteristic shapes are illustrated in Fig for an object with two plane 
faces. When the curvature is towards such a face, it is concave, and the plasma flows into the 
probe with a negative parallel velocity greater than 1 — My cot 6. This concavity is avoided 
in general if the probe has a convex curvature that is greater than 1/Ln- The characteristics 
curve either away from or towards the line x = const, depending on whether L„ is positive 
or negative. {My is taken always positive.) 

To summarize, then, the ion flux per unit perpendicular area from the plasma to a 
convex surface in the presence of combined (colinear) EAB and density-gradient diamagnetic 
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external drifts can be written precisely as eq (fT5|) : 

Ty = ncs = noo{zoo)cs exp[-l - (Myoo - Mh cot 6)], (38) 

but with 

Mh = ME + Mr^i (39) 

and 

riooizoo) = n^{z - 6z) = noo{z) exp{-6z/L„). (40) 

Because of its near even parity in My + 1, the Sz term can generally be ignored for the 
purposes of M± determination. 

5 Effects from presheath displacement 

The expressions obtained so far are for the flux leaving the plasma region where the drift 
expressions hold. Between that region and the probe surface lie the magnetic presheath, of 
order a Larmor radius thick, which we assume is small compared with the probe, and the 
Debye sheath, of order 4 Debye lengths thick, which we shall ignore altogether. The dynamics 
in the magnetic presheath are not ignorable. The flux to the probe itself is different from the 
flux from the plasma into the magnetic presheath when diamagnetic drifts are important. 

5.1 Magnetic presheath displacement calculation 

In the magnetic presheath, the electric field normal to the surface is strong enough that 
E±/B is of order the sound speed. The normal gradients of the resulting drift give rise to 
non-negligible convective derivative (i.e. polarization drift) terms. Indeed, it is those terms 
that permit the ion fluid trajectory to acquire sufficient cross-field velocity to satisfy the 
Bohm condition at the Debye-sheath edge. We assume on the magnetic presheath scale 
the probe surface can be approximated as planar and the gradients can be considered all 
to be normal to it. Let the direction of the normal (outward from the plasma) to the 
probe surface be k. Define sin a = — k.B/i? and the unit vector I = B A k/Bcosa in the 
direction perpendicular to B and k. Define the third unit vector by j = k A 1. See figure 
[6l Denote vector components in the respective directions by subscripts. Then the (normal) 
k-component of the momentum equation ([2]) is 

mw.Vvk = -ZeVk4> - (1 /n)VkP + Ze\. (B A k) (41) 

we eliminate the potential gradient using eq (fT6l) . and density gradient using the continuity 
equation ([T]) in the form nvk = const., so that V^lnn = — Vfclnffc, to obtain 

viQi cos a = {-cl/vk + Vk)VkVk ■ (42) 

We wish to calculate the total displacement, 6i in the l-direction, experienced by the ion 
fluid traversing the magnetic presheath. This is simply the time integral of eq fH21) : 

6i^i cos a = I {-cl/vk + Vk)VkVkdt = {1 - cl/vl)dvk = [vk + c^/vk] , (43) 
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Figure 6: (Color online). Illustration of the coordinates referred to the field and drift direc- 
tions (x, y, z) and the unit vectors referred to the probe surface directions (j, k, 1). 



(recognizing that Vkdvk/dxk = dvk/dt). The limits of the integral are Vk/cg = Tk/nCg = S aX 
the magnetic presheath outer edge (where is the normal flux density) and Vk/cs = 1 at 
its inner edge, entering the Debye sheath (whose thickness we ignore). Therefore 

6i = ps[2 -S- l/S]/cosa = -p.^r^^ . (44) 

b cos a 

Thus, for small incidence angle of the magnetic field {S small) a rather large displacement 
{6i ~ Ps/S) along the magnetic presheath occurs. This displacement has been noted and 



roughly estimated [18| in previous discussions of magnetic presheath structure. Its signifi- 
cance is that ions exit the magnetic presheath (and are collected by the probe) a tangential 
distance 6i from where they entered it. If there is a tangential gradient V^F^ of the normal 
flux-density F^ entering the magnetic presheath at position xi, then the fiux-density to the 
probe will be not Ffc(x/) but Ffc(x/ — 6i) = Ffc(xi) exp[— In(Ffc)] (choosing consistently 
with our In n assumption, uniformity of V In F^ when surface curvature is ignored). This 
alteration of the fiux-density is precisely the phenomenon that Cohen and Ryutov [19\ cal- 
culated in the small-a limit. In their Eulerian viewpoint, the alteration is attributable to 
divergence of the tangential flux in the magnetic presheath. From the present Lagrangean 
viewpoint, it arises from the integrated convective derivative. We can demonstrate this by 
evaluation. Identifying V In F^. = V In Uoo and using the definitions of S and 1 we find 

TOlnF, = -Ii^p;i.VlnF, = --[l^^^^±^(VlnnooAB).k 
Dcosa Dcos'^ac, ZeB-^ 
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[1 - sr 

S cos^ a Co 



Vne)-k 



b cos^ a 



M„ 



(45) 



Cohen and Ryutov's calculation assuming small S and a gave this expression with the 
geometric term [1 — S]'^ky/ cos^ a equal to 1, yielding a SiViTk equal to the difference between 
the ion and electron diamagnetic normal ffux densities. Now we evaluate the total ffux 
density to the probe using eq ( !38l) . which yields S = sin a, and recognizing that when the 
orientation of k perpendicular to B is defined by setting kz = sin j3, it immediately follows 
that kx = — cos [3 sin 6', ky = cos /5 cos 9 and sin a = cos P sin 6*. We obtain 



Tk{xi — 5i) = ncs sin a exp 





■ 1 — sin a ' 


{- 


. 1 + sin a. 



{Mr,i- Mne) cote 



The ion flux to the probe surface per unit perpendicular area is then 

f r 1 — sin ol 1 

Tiip = UooCs exp <^ -1 - Miioo + Mh cot 9 - — — ■ (M„i - 

L Ll + sma. 

The logarithm of ratio of the flux for positive and negative cos 9 is 



M„,)cote 



(46) 



(47) 



re 1 * 

ln(r+/r_) = - +7 " + cot 9 - 2^^^^|^(M„i - M„e) cot . (48) 



Ln 

Ignoring the 5z term, this becomes 

in(r+/r_) _ 

2 cot 6* 



sm a 



Me + 



Mfij2sina + M„e(l — sin a) 
1 + sin a 



(49) 



Thus, in addition to Me, the velocity combination measured by the Mach probe ratio is 
an interpolation between the ion and electron diamagnetic velocities (which of course have 
opposite sign), dependent upon the angle a between the probe face and the field. Normally 
in practice an intermediate value of 9 (not too small but not too close to 7r/2) must be used, 
which generally means an intermediate value of a. 



5.2 Probe curvature and finite facet size 

The displacement in the magnetic presheath also gives rise to a flux correction when the 
probe surface has curvature. The flux at any point on an electrode is characteristic of the 
flux into the magnetic presheath at a position —6i away. If the surface is curved, then 
the angles at that position will be different, giving rise to a different flux. But also, the 
displacement may have a divergence, which gives rise to flux enhancement even if the flux 
into the magnetic presheath were uniform. 

We note that the displacement strictly contains a component Sj, along the direction 
j = k A i, which for small a is mostly along B, and is of approximate magnitude ps/ sin a, 
similar to 6i . This displacement Sj can be calculated in the form of a closed integral expression 



by solving the magnetic presheath equations using the techniques of references |22l . |J]. We 
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ignore it, here and in the previous section, because it is mostly motion along the magnetic 
field, which can be considered to be accounted for by surface projection along the field, 
and because both the projection and any additional cross-field motion (in the k-direction) 
that makes 6j different from the pure field-line projection, make fiux contributions that are 
symmetric: of even parity under reversals of ky and hence of cos^. In other words, there 
is a small correction to the total flux density from 6j, which can be pictured as arising 
from the fact that the probe collects ions from a cross-fleld area that is larger than its solid 
cross-section by a margin of width approximately ps- But that correction does not affect the 
deduced Mach-numbers, because they are based on ratios of collection fluxes from surfaces 
with opposite ky, which are equally perturbed by even-parity terms. 

The perpendicular flux-density perturbation arising from 6i along 1 has contributions 
Sidn/dxi = SinMhd cot 9 / dxi and from the convective derivative of the perpendicular area 
element, which can be written (differentially) AA/A = V^{6il) = ^/V^.l + l.V*(5z, where V* 
denotes the two-dimensional gradient {d j dy,d / dz) in the perpendicular coordinates, but 
evaluated along the probe surface (not at constant x). Thus the total correction arising from 
surface curvature is 

- Ar||/r|| = Mh5i\.V' cot 9 + 5iV'\ + \.V'5i . (50) 

In this equation the flrst term has even parity with respect to fc^-reversal and therefore does 
not contribute to flow-measurement bias. The last two have odd parity and do contribute. 
Their order of magnitude is 5i/Rc where Rc is the typical radius of curvature of the sur- 
face. (Note that this is a real physical effect, not the elementary mathematical integration 
discussed by Pelemann et al|23|.) 

A related, but more intuitive and probably more relevant correction comes from edges 
of plane facets. Most Mach probes have electrodes at or near their ends (l^ . i H The 



direction 5J associated with the surfaces points forward, beyond the probe end, on one side 
and backward towards the probe stalk on the other, as illustrated in Fig [71 The surface 
for which it points backward collects essentially zero flux for a distance ~ 5i from the end, 
because of the magnetic presheath displacement; then it collects the full flux from then on. 
The surface for which 5 J points forward collects full flux at its end. The electrode will 
collect full flux throughout its area, provided that it is embedded in a surface of constant 
angle which extends a distance ~ 5i past the (backward) edge of the electrode. Therefore the 
likely effect (dependent on the detailed electrode placement) of 5i is to induce a depleted flux 
region on just the forward edge of just the electrode for which 6il is backward. An appropriate 
estimate of the odd-parity part of the resulting flux change is 

Arii/r,, ^ 6l/2h = _pA1^^1^ ^ , (51) 

2hS cos a 2h cos a sin 9 

approximating S as sin a to lowest order in M^, and denoting the average z-extent of the 
electrode as h. 

Since in a perpendicular Mach-probe measurement Mh is estimated from the expression 

Mh\ cot ^1 = (1/2) ln[r||+/r||_] , (52) 
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Figure 7: (Color online). Illustration of the opposite directions of magnetic presheath dis- 
placement {6il) for opposite values of cos^. The electric field is towards the probe and 5/1 is 
in the E A B direction. When the displacement is towards the stalk, as it is for facet 1 but 
not facet 2, a depleted region exists at the probe end. (In this figure "forward" is up and 
"backward" is down.) 



(where subscript ± refers to the sign of cos^) the perturbation AM^ to the deduced Mh is 
related to the odd-parity perturbation to flux, Arn, by 



AMh = |tan^|Ar||/r|| 



Ps (1 



sma 



2h cos a cos 9 



(53) 



There is therefore an intrinsic bias of order Ps/2h in such a Mach-probe measurement. Its 
direction is such that if the gradient of plasma pressure is in the forward direction of the 
probe, in other words the probe is introduced from the "outside" of the plasma where pressure 
is low (as is generally the case), then the spurious apparent drift is in the same direction as 
the ion-diamagnetic drift. This end-effect can in principle be avoided if the electrode does 
not sample the plasma at the end of the facets, but instead approaches the end no closer 
than a z-distance of approximately Ps{l — Sy cos j3 / {S cos a) . 



6 Temperature gradient drifts 

Diamagnetic drifts might also arise from temperature gradients, which have so far been 
excluded. A physically justifiable approach for electron temperature gradients simply regards 
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Te as constant along the field, but having an externally imposed gradient in the 2;-direction. 
A similar mathematical ansatz will be applied to Tj but with less clear physical justification. 
The equations of continuity and parallel momentum ([HD may once more be rearranged into 
characteristic form to obtain 



(v.V ± c^V)[ln?2/?2oo ± v\\/cs] = ±{v\\/cs){-vjLc) - vjLp 



(54) 



Here there are two inhomogeneous terms on the RHS. The first arises from the variation of 
Cg with z because of temperature gradients. Its gradient has been written as dcs/dz = Cs/Lc- 
The second arises because of (a) the gradient of Uoo (compare with eq (1371) ) and (b) the 
divergence v_l acquires when there are temperature gradients, which can be written V.v_l = 
Vz/Lt, where Lt = Ti/{dTi/dz) is the scale-length of ion-temperature gradient. Because (a) 
and (b) are of identical form, we combine them as —Vz/Lp with 1/Lp = 1/ Lt + 

The inhomogeneous terms on the RHS are the same order as the LHS terms near the 
boundary of the unperturbed region. They cannot therefore simply be ignored because 
of ordering. Their effects are to change the value of Inn/rioo — v\\/cs along the negative 
characteristic, which is what determines the fiux at the boundary. This change can be 
calculated, and proves to be usually ignorable, as we now show. We assume that the value 
of this negative characteristic quantity can be expressed as a function only of Mn: 



lnn/noo-M|| = g{M\\) 



(55) 



and solve for g using the same approach as we used for density gradients. In outline, the 
calculation is as follows. 

The positive characteristic equation gives 



+ 



Lr. 



We define dg/dM\\ = 2q so that ^ 
characteristics to obtain 

d /-Mil 

dt ^ = -"^ Ur 



d 
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(21nn/noo - g) 



(56) 



Mil and eliminate Mm derivatives from the 



± 



+ 



d_ 
di 



+ 



Lp l + q 



We eliminate partial x-derivatives between the positive and negative characteristics: 



- - 1) 



9 + (M|| + 1) 
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2csMhT-g 
dy 



(57) 



(5^ 



We re-express y-derivatives using the 2;- velocity expression fl27j) . which yields 

Mz = psil + l/2q)dg/dy. (59) 
Substituting (IFTI) and (IS^l) into (1351) and dividing the resultant through by Vz we arrive at 



- (^11 




l+q 



-M|| 
1^ 



+ 



p,(l + l/2g) 



(60) 
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This quadratic equation for q can be solved to obtain 



1 dg 



2(Mf 



s-w) + {l- s)M\\ ± ,/[(! - s)2 - 4]M|? + 2(s + 2)w 



2dM\\ 



-4(Mf 



w) 



(61) 



where w = MhLd Ps and s = Lc/ Lp. 

Provided M\\oo, w, and s are indeed constant, this is a consistent solution. In view of 
the variation of Cg and hence ps with z, it is clear that the constancy of u requires Cs{z) to 
satisfy a simple differential equation, whose solution gives the required Cs{z). This is not as 
simple a profile as the linear Innoo profile for the density-gradient. But we are interested in 
a local region on the gradient, where it is acceptable to choose the precise form of Cs{z) to 
satisfy the consistency condition without substantially changing the problem, in-so-far as a 
local expansion of Cg is allowable. 

The limit of eq (16T]) . when s and w — > oo together, is minus eq (!35l) although with different 
notation. Equation (ETil . can also be integrated analytically, but the resulting expression is 
long and cumbersome, so it is not given here. Instead, numerical integrations are illustrated 
in Fig [HI We choose to illustrate s = 0, which corresponds to a case where only the electron 
temperature has an external gradient. The case Mmoo = 0.2 is shown, to emphasize the 
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Figure 8: Solutions for the quantity g + M\\oo = Inn/noo — M\\ + M\\^ accounting solely for 
perpendicular gradients of sound-speed (s = 0). 



applicability of the solution shapes to all Mnoo- The solution becomes imaginary for some 
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relevant negative values of My when < w ^ 7, illustrated by the w = 1,3 curves. Provided 
such parameters are avoided, it can be seen that the values and especially the odd parity (in 
M|| + 1) part of g + Myoo are small. 

A second modification of the governing equations, arising from temperature gradients, is 
that additional terms are present in the drift velocity eq Q. Ion temperature gradients are 
easy to treat because they can be expressed as a uniform term like V0oo; but Tg gradients 
introduce a new type of term (c.f. [25]) . The total drift velocity is: 



V_L = - 



Tf.\ „ / , Ti \ Tj mc 



ln(n/noo) V ( 7^ j + ^ + ^ j + In "^^o + ^ V In n/n. 



2 



A^- (62) 



We can still discard the final (Vlnn/noo) term as an M^i giving no advection because 
M^i.Vlnn/noo = 0. The rest should be regarded as Vh, having no 2;-component. The first 
term is the new one, so far not treated. It is non-uniform; but a function of Inn/uoo and 
hence, within the solution schemes we have developed, a function of My. The Vh variation 
can be incorporated into the scheme above, except that the consistency condition upon the 
Cs profile becomes that w (= MhLc/ps), which varies with Vh, must be a function only of My. 
Assuming an appropriate Cg shape is chosen, eq ( I6T1) still applies, and could still be integrated 
provided that w{M\\) were known. Plainly, provided that w stays within acceptable ranges 
during the integration, curves that bear resemblance to those of Fig ([8]) will be obtained. 
Generally then, the arguments that g + M\\oo can usually be ignored will remain as valid as 
before, provided w stays within the range of their validity. 

The boundary condition at the plasma edge is still that the positive characteristic be 
tangent to the probe surface, that is (eq [14]) M\\ = Mh cot 9 — 1, but with the following 
definition: 



ln(n/noo)V (^^) + ^ (|^ + '^~) + ^^^^'^'^ 



- \n{n/n^)MTe + Mn + Me + M^i , (64) 



where Mxe and M^j are the (external) diamagnetic drift due to temperature gradient of the 
electrons and ions. We can then substitute into eq fl55l) to eliminate My and obtain after 
rearrangement: 

-1 + g + {MTi + Me + M^i) cote 
mn Tioo = , (do) 

' I +MTe cote ' ^ ' 

in which, for the purposes of perpendicular velocity measurement, we can take g ~ —M\\^. 
To first order in the perpendicular velocities this can be written 

Inn/n^ ^ -1 - M\\^ + [(1 + M\\^)MTe + + Me + M^i] cot 9 , (66) 

showing that the effect of the electron temperature-gradient drift can be amplified or atten- 
uated, compared with the other drifts, depending upon the parallel external velocity Myoo- 
Finally we must account for the temperature gradients in the correction arising from 
magnetic presheath displacement. We can ignore the next order corrections to 5i arising 
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from approximations in its derivation. But we must account for the fact that since 6i oc c^, 
a transverse derivative of Cg gives rise to transverse divergence which alters the flux density 
as it traverses the magnetic presheath. The change in area A arising from this effect is 
AA/A = SiLVlnCg. This gives rise to a change in InF of — ^J.VlnCs. In addition, the 
displacement gives rise to a convective difference between the flux to the probe and that 
entering the presheath that, since F oc ncg can be expressed as 



A\n{ncs) = -6il.{Vlnn + Vlnc^) 
adding these two effects gives the total magnetic presheath difference: 

AlnTfc = -5ri-(Vlnn + 2Vlnc,) = -6il.(y Inn + V \n{ZTe + Ti)) 



(67) 



(68) 



Applying the transformations of eq (145|) . we see that our result is the same as before except 
that the total diamagnetic difference velocity Md = M^i — MDe = {Mni + Mti) — {Mne+ Mxe) 
is involved rather than just the density- gradient part calculated before. 

The final expression for flux per unit perpendicular area to the probe is thus 
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A + g + {Mn + Me + M^i) cot 6 
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sma 
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sma 



1 + Mtc cot 9 

or, to first order in perpendicular velocities, and approximating g, 

' 1 — sin a 



{Moi - Moe) cote), (69) 



-1 - Aflloo + 



;i + M||oo)MTe + Mm + Me 



1 + sin a 



D 



cot 9, (70) 



where a is the angle between the probe surface and the magnetic field (in 3-dimensions), 
and 9 is the angle within the plane containing field and external drift. All drift velocities 
here refer to the external, unperturbed plasma. 



7 Discussion 

The complete solution of the drift equations that has been obtained here is highly appropriate 
when the E A B drift dominates. Then the positive characteristics are straight and the full 
solution in the plasma region can readily be constructed. Notice that the parallel length of 
the presheath (i.e. the perturbed plasma region) is approximately the transverse size of the 
probe Tp (say) times (My + 1)/M±. If the perpendicular Mach number is not very small, 
this length ~ rp/M± is likely substantially shorter than that from the standard diffusive 
estimate: r^Cs/D. Certainly it can easily be shorter than the mean-free-path for electron-ion 
Coulomb collisions, which is required for the ignoring of friction inherent in our treatment. In 
such a situation taking the electron temperature to be invariant along the field is completely 
natural. It is not so natural to make that approximation for the ions. However, it is 



known from other calculations [25|, l9| that the isothermal approximation gives results quite 



close to those that arise from more physically plausible approximations. In any case, the 
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standard widely-accepted formulas are based upon isothermal-ion calculations. In addition to 
providing rigorous analytical justification, regardless of probe geometry, for formulas that are 
practically the same as those arising from a diffusive treatment, the present treatment helps 
to resolve another conceptual problem of long standing. It is that probes are often smaller 
than the typical transverse size of the turbulence that is responsible for transport in the 
regions of plasma in which they are used. In other words, transport, for example in tokamak 
edges and scrape-off-layers, is actually known to be dominated by fluctuating cross-field 
flows that in many situations have eddies larger than the probes. In such situations, using a 
heuristic approximation that cross- field flux is expressible through a diffusion coefficient, as 
prior treatments have done, is questionable. The present treatment, regarded as a short-time 
snap-shot of a situation that is fluctuating, is more appropriate. Fortunately the result is 
the same, although with much sharper physical justification. 

The drift-approximation is unproblematic for E A B drifts, because there is nothing to 
prevent the probe being bigger than the Larmor radius while the drift Mach number is of 
order unity. This is not so for diamagnetic drifts. The diamagnetic Mach number is of 
order Md ps/L (where L is the pressure scale- length). Consequently if Md ~ 1, there 
is no separation of scales between the Larmor radius and the gradient scale-length. It is 
then impossible to choose a probe size that is both bigger than the Larmor radius (in order 
to justify the drift approximation) and smaller than the gradient scale-length (to justify a 
local approximation of the drifts). This is an inherent difficulty that underlies the need 
in the calculations to choose a specific shape of the plasma profiles. Only if Mo is small 
does a local approximation to the flow make sense. Therefore, the usefulness of the present 
calculation is increasingly compromised as Mo approaches unity. A calculation that avoids 
the drift-approximation is then really needed. It seems unlikely that an analytic solution 
like the present one will be forthcoming. 

In summary, complete solutions of the problem of ion collection by an arbitrary-shaped 
object have been obtained in the drift approximation (ignoring the polarization drift). The 
normalized flux density fIBIJl) . (170]) is a function only of the orientation of the surface, provided 
the object is convex. The precise meaning of "convex" here is that no positive characteristic 
that originates elsewhere on the object should pass through the point of interest. A transverse 
Mach probe using a variety of electrode orientations measures the E A B drift plus this 
specified combination of both ion and electron diamagnetic drifts. Its ideal calibration has 
been derived, and possible problems arising from finite size identified and quantified. 
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